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We study the existence and properties of rogue wave solutions in different nonlinear wave evolution 
models that are commonly used in optics and hydrodynamics. In particular, we consider Fokas- 
Lenells equation, the defocusing vector nolinear Schrodinger equation, and the long-wave-short- 
wave resonance equation. We show that rogue wave solutions in all of these models exist in the 
subset of parameters where modulation instability is present, if and only if the unstable sideband 
spectrum also contains cw or zero-frequency perturbations as a limiting case (baseband instability). 

We numerically confirm that rogue waves may only be excited from a weakly perturbed cw whenever 
the baseband instability is present. Conversely, modulation instability leads to nonlinear periodic 
oscillations. 
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I. INTRODUCTION 

Many nonlinear wave equations associated with differ¬ 
ent physical systems exhibit the emergence of extreme, 
high-amplitude events that occur with low probability, 
and yet may have dramatic consequences. 

Perhaps the most widely known examples of such pro¬ 
cesses are the giant oceanic rogue waves [1] that unex¬ 
pectedly grow with a great destructive power from the 
average sea level fluctuations. This makes the study of 
rogue waves a very important problem for ocean liners 
and hydrotechnic constructions mm- Hence, it is not 
surprising that the phenomenon of rogue waves has at¬ 
tracted the ample attention of oceanographers over the 
last decade. Although the existence of rogue waves has 
been confirmed by multiple observations, uncertainty still 
remains on their fundamental origins [4]. 

In recent years, research on oceanic rogue waves has 
also drawn the interest of researchers in many other do¬ 
mains of physics and enginering applications, which share 
similar complexity features: in particular, consider non¬ 
linear optics [5]. The ongoing debate on the origin and 
definition of rogue waves has stimulated the comparison 
of their predictions and observations in hydrodynamics 
and optics, since analogous dynamics can be identified 
on the basis of their common mathematical models [6]. 

So far, the focusing nonlinear Schrodinger equation 
(NLSE) has played a pivotal role as universal model for 
rogue wave solutions, boh in optics and in hydrodynam¬ 
ics. For example, the Peregrine soliton, first predicted 
as far as 30 years ago [7], is the simplest rogue-wave so¬ 
lution of the focusing NLSE. This rogue wave has only 
recently been experimentally observed in optical fibers 
[8] , water-wave tanks [9] , and plasmas [10] . 

Eor several systems the standard focusing NLSE turns 
out to be an oversimplified description: this fact pushes 
the research to move beyond this model. In this direction, 
recent developments consist in including the effect of dis¬ 


sipative terms. In fact, a substantial supply of energy 
(f.i., from the wind in oceanography, or from a pump¬ 
ing source in laser cavities) is generally required to drive 
rogue wave formation m- Because of their high ampli¬ 
tude or great steepness, rogue wave generation may be 
strongly affected by higher-order perturbations, such as 
those described by the Hirota equation [12], the Sasa- 
Satsuma equation m and the derivative NLSE m- 

The study of rogue wave solutions to coupled wave sys¬ 
tems is another hot topic, where several advances were 
recently reported. Indeed, numerous physical phenom¬ 
ena require modeling waves with two or more compo¬ 
nents. When compared to scalar dynamical systems, vec¬ 
tor systems may allow for energy transfer between their 
different degrees of freedom, which potentially yields rich 
and significant new families of vector rogue-wave solu¬ 
tions. Rogue-wave families have been recently found as 
solutions of the vector NLSE (VNLSE) [T5HT8] . the three- 
wave resonant interaction equations [19] , the coupled Hi¬ 
rota equations [20], and the long-wave-short-wave reso¬ 
nance [21]. 

As far as rogue waves excitation is concerned, it is 
generally recognized that modulation instability (MI) 
is among the several mechanisms which may lead to 
rogue wave excitation. MI is a fundamental property 
of many nonlinear dispersive systems, that is associated 
with the growth of periodic perturbations on an unstable 
continuous-wave background [22] . In the initial evolution 
of MI, sidebands within the instability spectrum experi¬ 
ence an exponential amplification at the expense of the 
pump. The subsequent wave dynamics is more complex 
and it involves a cyclic energy exchange between multi¬ 
ple spectral modes. In fiber optics, MI seeded from noise 
results in a series of high-contrast peaks of random in¬ 
tensity. These localized peaks have been compared with 
similar structures that are also seen in studies of ocean 
rogue waves [5]. Nevertheless, the conditions under which 
MI may produce an extreme wave event are not fully un- 
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derstood. A rogue wave may be the result of MI, but con¬ 
versely not every kind of MI necessarily leads to rogue- 
wave generation [181 l23ll25] . 

In this work, our aim is to show that the condition for 
the existence of rogue wave solutions in different nonlin¬ 
ear wave models, which are commonly used both in optics 
and hydrodynamics, coincides with the condition of base¬ 
band ML We define baseband MI as the condition 
where a cw background is unstable with respect 
to perturbations having infinitesimally small fre¬ 
quencies. Conversely, we define passband MI the 
situation where the perturbation experiences gain 
in a spectral region not including cj = 0 as a limit¬ 
ing case. We shall consider here the Fokas-Lenells equa¬ 
tion (FLE) [14], the defocusing VNLSE [18] and the long¬ 
wave-short-wave (LWSW) resonance [21]. As we shall 
see, in the baseband-MI regime multiple rogue waves can 
be excited. Conversely, in the pass-band regime, MI only 
leads to the birth of nonlinear oscillations. 

We point out that, in this work, we consider 
as rogue wave a wave that appears from nowhere 
and disappears without a trace. More precisely, 
we take as a formal mathematical description of 
a rogue wave a solution that can be written in 
terms of rational functions, with the property of 
being localized in both coordinates. 


II. FOKAS-LENELLS EQUATION 

The FLE is partial differential equation that has been 
derived as a generalization of the NLSE IMIEU- In the 
context of optics, the ELE models the propagation of 
ultra-short nonlinear light pulses in monomode optical 
fibers [27] . 

Eor our studies, we write the ELE in a normalized form 


i(l -h indr)^^ + + cr|'0p(l + iKdr)^ = 0, (1) 


where represents the complex envelope of the 

field; r are the propagation distance and the retarded 
time, respectively; each subscripted variable in Eq. 0 
stands for partial differentiation, a {a = ±1) denotes a 
self-focusing {a = 1) or self-defocusing (a = —1) non¬ 
linearity, respectively. The real positive parameter k 
{k > 0) represents a spatio-temporal perturbation. Eor 
K, = 0, Eq. Q reduces to the NLSE. 

Soliton, multi-solitons, breathers and rogue waves so¬ 
lutions have been recently found for Eq. Q. Let us 
examine the existence condition for these rogue waves. 
The rogue wave solutions may be expressed as na 




2iK^{^ + 2kt) + aK/a^) 

D + ihiK^ 


( 2 ) 


where 'i/jq = represents the background solu¬ 

tion of Eq. 0> a is the real amplitude parameter (a > 0), 


uj the frequency; moreover [3 = j2K —go?^ K = 

gK/{AP), a = ±VgK - 


The rogue wave solutions 0 depend on the real pa¬ 
rameters a and cj, for fixed a and k. In the focus¬ 
ing regime (a = 1), rational rogue waves exist for uj 
in the range [1//^ — 1/{Pk^)P/k]. Whereas in the de¬ 
focusing regime rogue waves exist for uj in the range 
[I/kP/k + Eigure shows the domains of 

rogue wave existence in the plane {uj,k)^ for either the 
focusing or the defocusing regimes. Surprisingly, expo¬ 
nential soliton states exist in the complementary region 
of the (cj, k) plane (see Ref. la for details on the prop¬ 
erties of these nonlinear waves). Eigur e illustrates a 
typical example of rogue wave solution ([2fT^ 



FIG. 1. Existence domains of rogue waves in the plane (ac, uj), 
with a = 1, in the focusing regime (cr = 1) and defocusing 
regime (cr = — 1). The red dotted line denotes uj = 
green dashed line u = 1//^; black solid line u = 1/ac — 1/P. 


Let us turn our attention now to the linear stability 
analysis of the background solution of Eq.Q. A per¬ 
turbed nonlinear background can be written as pp = 
[a where p(^, r) is a small complex pertur¬ 

bation that satisfies a linear differential evolution equa¬ 
tion. Whenever p is r-periodic with frequency Q, i.e., 

such equation reduces 
to a set of 2 X 2 linear ordinary differential equations 
rj' = iMf], with rj = (here a prime stands for 

differentiation with respect to r). Eor any given real fre¬ 
quency Q, the generic perturbation r]{^) is a linear com¬ 
bination of exponentials where Wj,{j = 1,2) are 

the two eigenvalues of the matrix M = whose 

elements read as: 
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G(Q) = —Im{iL’} > 0 , perturbations grow exponentially 
like exp(G^) at the expense of the pump wave. 

MI is well depicted by displaying the gain G(fl) as func¬ 
tion of n and Q. The resulting MI gain spectrum 

is illustrated in Fig. [^and Fig. These figures show 



FIG. 2. Typical rogue soliton states. Top, focusing regime 
(7 = 1, = 0.5 and a = l,ct; = 0. Bottom, defocusing regime 

a = —1, K. = 0.5 and a = l,ct; = 4. 





Mil 

Mi 2 
M 21 
M 2 2 


— + (jo^K — VL{uj + + (jo^n) 

{K - kQ) 

(jo?K 
{K - kQ) ’ 
aa^K 

— (to^K — Vt(uj + + (70? n) 

{K + kQ) ' 


Since the entries of the matrix M are all real, the eigen¬ 
values Wj are either real or they appear as complex con¬ 
jugate pairs. The eigenvalues of the matrix M are the 
roots of its characteristic polynomial, 

B[w) = B 2 W^ + Biw + Bq^ (3) 

B 2 = 

Bi = —4Q{2I3 kK + 2Kuj)^ 

Bq = 

2Pi<i{o? tier + cj) + A^a{K + 2k,uj). 

Mi occurs whenever M has an eigenvalue w with a 
negative imaginary part. Indeed, if the explosive rate is 


FIG. 3. Maps of logaritmic MI gain (lOlogioG) in the focusing 
(cr = 1) FLE 0. Top, MI on the (Q, o) plane, calculated for 
the case a = — 1. Bottom, MI on the (Q,ct;) plane, 

calculated for the case a = 1, /c = 0.5. 

the MI gain in the focusing ans defocusing regime, respec¬ 
tively. In both cases, baseband MI is only present in a 
certain subset of the cj, n parameters. Since the gain band 
(where G{Q) 7 ^ 0 ) can be written as 0 < Qi < Q < Q 2 
(and its symmetric counterpart with respect to = 0 ), 
baseband MI is obtained if f2i = 0, whereas passband MI 
occurs for fti > 0 . 

We proceed next by focusing our attention on the MI 
gain spectrum, by evaluating the sign of the discriminant 
A of the characteristic polymomial (§: this leads to 

sign{A} = sign{ft^ — Aa^aK^ (l — (4) 

If the discriminant A is positive, the characteristic 
polynomial has two real roots and there is no MI. On 
the other hand if the discriminant A is negative, the 
characteristic polynomial B has two complex conjugate 
roots, and Eq. 0 exhibits baseband MI. It is clear from 
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FIG. 4. Maps of of logaritmic MI gain (10/opioG) in the defo- 
cusing (cr = —1) FLE Q. MI on the {Q,k,) plane, calculated 
for the case a = 2, ct; = 1. 


Eq. Q that for FLE if there is MI, it is of baseband 
type only: either the system is modulationally unstable 
for O ^ 0, either there is no MI at all. The interest¬ 
ing finding is that the sign constraint on the discrim¬ 
inant, which determines the presence of baseband MI, 
leads to the condition that uj should be in the range 
[1/k — l/{a^ k] in the focusing regime (a = 1 ), 
and in the range [ 1 //^, l/zi: + l/(a^/^^)] in the defocusing 
regime (a = — 1). These conditions exactly coincide with 
the constraints that are required for the existence of the 
rogue wave solution <§■ 

These results are important since they show that, for 
both the focusing and the defocusing regime, rogue wave 
solutions of Eq. Q only exist in the subset of the pa¬ 
rameters space where also baseband MI is present. 

We checked the results of our analysis by extensive 
numerical solutions of Eq. 0 - These simulations in¬ 
deed confirm that, in the baseband MI regime, multi¬ 
ple rogue waves can generated from an input plane wave 
background with a superimposed random noise seed (see 
Fig. §. 


III. DEFOCUSING VNLSE 

The defocusing VNLSE constitutes another model that 
has been thoroughly exploited for the description of fun¬ 
damental physical phenomena in several different disci¬ 
plines. In oceanography, for instance, it may describe 
the interaction of crossing currents [28]. In the context 
of nonlinear optics, it has been derived for the description 
of pulse propagation in randomly birefringent fibers [29| , 
or coupled beam propagation in photorefractive media 
[SO]. 

For our studies, we write the defocusing VNLSE in the 
following dimensionless form 


EIG. 5. Golor plot of t)|^ from the numerical solution of 
the focusing ELE 0 in the baseband MI regime. The initial 
condition is a plane wave perturbed by a random noise seed, 
with parameters: a = 2, cj = 1, = 1.15, a = —1. 


-h -h = 0 

-h P -h = 0, 

where represent complex wave en¬ 

velopes; ^,r are the propagation distance and the re¬ 
tarded time, respectively; each subscripted variable in 
Eqs. 0 stands for partial differentiation. Note that Eqs. 
0 refer to the defocusing (or normal dispersion) regime. 
Unlike the case of the scalar NLSE, rational rogue solu¬ 
tions of the defocusing VNLSE do exist, as it was recently 
demonstrated [18]. These rogue wave solutions can be 
expressed as: 

J,(j) = ,1,0) + pTjaj +130j) - + pOj . 

^ ^ + ;5(pr + 1 ) 

(6) 

with j = 1 , 2 . represent the back¬ 

ground solution of Eqs. 0, aj are the real ampli¬ 
tude parameters {aj > 0 ), uoj are the frequencies, and 

Moreover, aj = 4 p^/(p^ + ~ ~ 

ip); 13 = p^/x{p‘^ + 4 cjiu;2 ),P = 2 Im(A + k),uji + u;2 = 
2 Re(A k) ^ uji — UJ2 = 2ci;, y = Imk. The evaluation of the 
complex value of A and k should be performed as follows. 
The parameter A is the double solution of the polynomial 
M(A) = A^ -|- M2A^ T M]^A “I" Mq = 0 , with Aq = — k^ -|- 
k{uj‘^-\-ai-\-a 2 )-\-uj{a 2 — ai), Ai = —cj^-hai+al, A2 = 
k. Moreover, the constraint on the double roots of A{X) 
is satisfied whenever the discriminant of ^(A) is zero, 
which results in the fourth order polinomial condition 
D{k) = k^ Dsk^ + D2k‘^ + Dik -f- Do = 0 , with Dq = 
(cj^ —a\ — 02 )^/ (2^co’^) — (3/4)^(a2 — Di = —9(a2 — 
a\){2uj‘^ a\a 2 )/{ 2 ^uj)^ D 2 = ~[8^^ — {a\ + <^ 2 )^ T 
20 co’^(a^ + a2)]/(2^u;^), D3 = (a^ — a^)/(2ci;). Thus, A is 
the double solution of the third order polynomial A{X), 
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and k is any strictly complex solution of the fourth order 
polynomial D{k) (see Ref. [18] for details on nonlinear 
waves calculations and characteristics). 

The rogue waves (§ depend on the real parameters 
ai , a 2 and uj which originate from the backgrounds: ai , a 2 
represent the amplitudes, and 2 uj the “frequency” differ¬ 
ence of the waves. Figure shows a typical dark-bright 
solution 



FIG. 6. Rogue wave envelope distributions ^)| and 

(f)| of expression Here, ai = 3,a2 = 1,0; = 1. 
k = 2.36954 -h 1.1972z and A = -1.69162 - 1.79721z. 

In the defocusing regime, it has been demonstrated 
HHI that rogue waves exist in the subset of parameters 
ai, a 2 , CO where 

(a^+a^)^ —12(a^—7a^a2+a2)co^+48(a^+a2)co^—64co^ > 0. 

(7) 

Figure illustrates two characteristic examples of the 
existence condition for rogue waves. In particular, Fig0 
shows that, for a fixed co, the background amplitudes 
should be sufficiently large in order to allow for rogue 
wave formation. 

Let us turn our attention now to the linear stability 
analysis of the background solution of Eqs.([^. A per¬ 
turbed nonlinear background may be written as 'ipp^ = 
[aj where are small complex per¬ 

turbations that obey a linear partial differential equa¬ 
tion. Whenever Pj{^,r) are r—periodic with frequency 




FIG. 7. Rogue wave existence condition, a) (cu^ai) plane, 
with a 2 = 3. b) (a 2 ,ai) plane, with cj = 4. 

Q,, i.e., Pji^,T) = , their equa- 

tions reduce to the 4x4 linear ordinary differential equa- 
tion 7]' = iMrj, with ij = ■ For 

any given real frequency the generic perturbation r]{^) 
may be expressed by a linear combination of exponentials 
ex.p{iwj^) where Wj^ j = !,••• ,4, are the four eigenval¬ 
ues of the matrix M = {Mij}. 

Mil = - 2Quji - 2 ai, 

M 22 — — 2f}c(Ji T 2(2^, 

M 33 = - 2^002 - 2al, 

AI 4 . 4 . = — 2^U)2 T 2(22: 

Mi2 = —M 21 = —2a^, 

Mi 3 = Mi 4 = M31 = M32 = —2aia2, 

M 41 = M 23 = M 24 = M 42 = 2 aia 2 , 

A743 = —Af34 = 2(22- 

Since the entries of the matrix M are all real, the eigen¬ 
values Wj are either real or they appear as complex con¬ 
jugate pairs. These eigenvalues are the roots of the char¬ 
acteristic polynomial B{w) of the matrix M: 

B{w) = Bsw^ + B 2 W^ Biw + 5o, 

Bo = - 4cj2)[4(a? + a" - 

Bi = 16uj{ai — 

B 2 = -2[2{al + + 2 cj^) + 

^3 = 0. 

MI occurs whenever M has an eigenvalue w with a nega¬ 
tive imaginary part, Im{u;} < 0. Indeed, if the explosive 
rate is G{Q) = —Im{u;} > 0, initial perturbations grow 
exponentially as exp(G^) at the expense of the pump 
waves. Typical shapes of the MI gain G{Q) are shown in 
Fig. [ 8 ] 

Figure |^a) corresponds to the case where the nonlin- 
ear background modes have opposite frequencies (cji = 
—UJ 2 = ^)- The higher cj, the higher G. In the special 
case of equal background amplitudes ai = a 2 = a, the 
marginal stability conditions can be analytically found: 

= maxldcj^ — 8 a^, 0}. Thus, for j2 
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Thus we may conclude that in the defocusing regime, 
rogue wave solutions (§ only exist in the subset of the 
parameter space where MI is present, and in particular 





FIG. 8. Maps of MI gain 2G of the VNLSE a) MI on 
the (Q,ct;) plane, calculated for the case ai = 3, a 2 = 3, 
(jUi = —(jU 2 = Ct;. Dotted (green online) curves represent 
the analytical marginal stability condition Q. — 2(jj^ 0?“ — 
max{4ct;^ — 8a^, 0}. b) MI on the (Q, ai) plane, calculated for 
the case a 2 = 3, ct;i = —002 = 4. 


a baseband MI, which includes frequencies that are ar¬ 
bitrarily close to zero, is present (i.e. 0 < < 4a;^). 

Instead, for < Ci;^/2, MI only occurs for frequencies 
within the passband range — 8 a^ < < 4u;^. We 

may point out that the rogue waves necessarily exist 
for > uj‘^/2. Thus, rogue waves (| 6 ) and baseband MI 
coexist. 

Figure [^b) illustrates the case of different frequencies 
{uji = —co ’2 = cj) and input amplitudes ai 7 ^ a 2 for the 
nonlinear background modes. For low values of ai, only 
passband MI is present. By increasing ai, the baseband 
MI condition is eventually attained. 

In order to analytically represent the condition for the oc¬ 
currence of baseband MI, let us consider the limit (2^0. 
To this aim, we may rewrite the characteristic polyno¬ 
mial as B{Qv) = and consider the polynomial 

b{v) at = 0 , namely b{v) = b^v^ + 62 ^^ + biv + 60 , 
bo = + a2 — bi = 16 g(a^ — a^), 62 = 

—4(ai + al + 63 = 0. Let us evaluate now the 

discriminant of the characteristic polynomial B: if the 
discriminant is positive, B has four real roots, and no 
MI occurs. Whereas if the discriminant of B is nega¬ 
tive, there are two real roots and two complex conjugate 
roots, and Eqs.([^ exhibits baseband MI. Again, the in¬ 
teresting finding is that the constraint on the sign of the 
discriminant of the characteristic polynomial 5, which 
leads to the baseband MI condition, turns out to exactly 
coincide with the sign constraint 0 that is required for 
rogue wave existence. 



EIG. 9. Golor plot of (a) and ^)| (b) from 

the numerical solution of the defocusing VNLSE. The initial 
condition is a plane wave perturbed by weak random noise. 
Parameters: ai = 2, a 2 = l,Cc;=l. A rogue wave is high¬ 
lighted by a surrounding box. 

Fig. |9] and Fig. show two different numerically 
computed nonlinear evolutions, obtained in the case of 
baseband MI (leading to rogue wave generation) and of 
passband MI, respectively. These evolutions permit to 
highlight that the nonlinear evolution of baseband MI 
leads to rogue wave solutions of the VNLSE ^ . Fig¬ 
ure shows the numerically computed evolution of a 
plane wave perturbed by a small random noise in the 
baseband MI regime. After a first initial stage of lin¬ 
ear growth of the unstable frequency modes, for ^ > 5 
the nonlinear stage of MI is reached. As we can see, MI 
leads to the formation of multiple isolated peaks (dips) 
that emerge at random positions. By carefully analyzing 
one of these peaks, for example the peak near the point 
(r = 0,^ = 9), we may clearly recognize the shape of a 
rogue wave as it is described by the expression Con¬ 
versely, FigfTol shows the numerically computed evolution 
of a plane wave perturbed by a small random noise, in 
the passband MI regime. After a first initial stage of lin¬ 
ear growth of the unstable frequency modes, for ^ > 2 
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FIG. 10. Color plot of <f)| (a) and ^)| (b) from 

the numerical solution of the defocusing VNLSE. The initial 
condition is a plane wave perturbed by weak random noise. 
Parameters: ai = 2,a2 = l,Cc; = 3. No rogue waves are 
generated in this case. 


where represents the short wave complex en¬ 
velope, and r) represents the long wave real field; 

^ and r are the propagation distance and the retarded 
time, respectively; each subscripted variable stands for 
partial differentiation. 




the nonlinear stage of MI is reached. In this case, we may 
observe the generation of a train of nonlinear oscillations, 
with wave-numbers corresponding to the peak of MI gain 
{^max = 5). As it was expected, no isolated peaks (dips) 
emerge from noise in this case, given that the condition 
for the existence of rogue waves is not verified. 


IV. LWSW MODEL 


FIG. 11. Rogue wave envelope distributions <f)| and 

(f)| corresponding to expressions Here, a = l,ct; = 
0,6 = 0.5. 

The fundamental rogue wave solution of Eqs. © has 
recently been reported in Ref. [21], and reads as 

j_ iz _^_1_ 

0 L ^ ^2^2 ^ ’ 


The last model we consider in our survey is the LWSW 
resonance. It is as well a general model that describes the 
interaction between a rapidly varying wave and a quasi 
continuous one. In optics the LWSW resonance rules 
wave propagation in negative index media m or the 
optical-microvave interactions [32| . Whereas in hydrody¬ 
namics the LWSW resonance results from the interaction 
between capillary and gravity waves [33] . 

For our studies, we write the LWSW equations in the 
dimensionless form 


-h = 0 


n^e-{r-m0^ + l/4n^ 

[(r - + n 2^2 + i/4n^]2 ’ 

where represents the background so- 

lution of the short wave, defined by the amplitude a 
{a > 0), frequency cj, and wave number P = l2 — h] the 
amplitude 6 (6 > 0) defines the background solution of 
the coupled long wave real field. The p arameters m an d 
n are real, defined by m = |[5ci; — y/3(cJ^'^rTH^^u7/)], 
n = (3m — uj)(m — cj), with v = + bcja^, 

P = I = -{p - yp2 _ y3)i/3^ 

for uj < —3(2a^)^/^, and I = {—p + p'^ — for 



( 8 ) 



















—3(2a^)^/^ < uj < |( 2 a^)^/^. LWSW rogue waves (|9 
depend on the real parameters a, uj and b (see Ref. 
for details on nonlinear wave characteristics). Figure 
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shows a typical LWSW rogue solution. Importantly, the 
existence condition for rogue waves of the LWSW model 
is that UJ < |( 2 a^)^/^. 


Let us turn our attention now to the linear stabil¬ 
ity analysis of the background solution of Eqs. 

Here a perturbed nonlinear background can be written 
as = [a ^ and ^ + Pl where 

are small complex perturbations that 
obey linear partial differential equations. Whenever the 
perturbations ps^PL are r—periodic with frequency 
i.e., ps{^,^) = recalling that 

is real, Pl{^^^) = 

turbation equations reduce to a 3 x 3 linear ordinary 
differential equation 77 ' = iMr]^ with 77 = 

(here a prime stands for differentiation with respect to 
r). For any given real frequency H, the generic perturba¬ 
tion may be expressed as a linear combination of expo¬ 
nentials ex.-p{iwj^) where tuj, j = 1, • • • ,3, are the three 
eigenvalues of the matrix: 


M = 


-^n^-ujn 

0 

fta 


0 

-ujft 

Vta 



( 10 ) 


Since the entries of the matrix M are all real, the eigen¬ 
values Wj are either real, or they appear as complex con¬ 
jugate pairs. These eigenvalues are obtained as the roots 
of the characteristic polynomial B{w) of the matrix M: 


B{w) = Bsw^ + B 2 w‘^ + Biw + Ho, (11) 

Ho = a^^l^ Hi = uj^n^ - 11^/4, H 2 = 2ujn, H 3 = 1 . 

MI occurs whenever M has an eigenvalue w with a 
negative imaginary part, i.e., Im{7n} < 0. Indeed, if the 
explosive rate is G(0) = —Im{7n} > 0, perturbations 
grow larger exponentially like exp(G^) at the expense of 
the pump waves. By calculating the discriminant of the 
polynomial H, one finds A = — uj{9a‘^ — 

uj^)Vt^ ^Ao?uj^ — 21 a^). If the discriminant A is positive, 
the polynomial H has real roots, and no MI occurs. Con¬ 
versely if the discriminant A is negative, the polynomial 
H has two complex conjugate roots, which means that 
MI is present for Eqs.(|^. The marginal stability curves, 
corresponding to A = 0, can thus be calculated. Eigure 
shows a typical MI gain spectrum of the LWSW Eqs. 
as one can see, there exist regions of either baseband 
or passband ML 

As in previous sections, let us proceed now to discuss 
the MI behavior in the limit situation where —> 0, a 
condition which characterizes the occurrence of baseband 
ML In this regime, the discriminant of the polynomial H 
reduces to A = Ao?uj^ — 27a^, which leads to the MI con¬ 
dition UJ < |(2a^)^/^. Again, the baseband MI condition 



FIG. 12 . Maps of MI gain 2G of the LWSW Eqs. MI on 
the (11, Ct;) plane, calculated for the case a = 1, Dashed (green 
online) curves represent the analytical marginal stability con¬ 
dition 11 ®(Tq 6 _ 1 q 4 _ _ 27 ^ 2 ^ _ 


turns out to exactly coincide with the condition for the 
existence of rogue wave solutions of Eqs. <§■ 

Eigure shows a numerical solution of LWSW, ob¬ 
tained in the case of baseband MI (leading to rogue wave 
generation), showing the evolution of a plane wave per¬ 
turbed by a small random noise. After a first initial stage 
of linear growth of the unstable frequency modes, for 
^ > 8 the nonlinear stage of MI is reached. As we can 
see, MI leads to the formation of multiple isolated peaks 
that emerge at random positions. By carefully analyzing 
one of these peaks, we may clearly recognize the shape 
of a rogue wave as it is described by the expression 


V. CONCLUSIONS 

In this work we studied the existence and the proper¬ 
ties of rogue wave solutions in different integrable nonlin¬ 
ear wave evolution models which are of widespread use 
both in optics and in hydrodynamics. Namely, we con¬ 
sidered the Eokas-Lenells equation, the defocusing vector 
nolinear Schrodinger equation and the long-wave-short¬ 
wave resonance. We found out that in all of these 
models rogue waves, which can be modeled as 
rational solutions, only exist in the subset of pa¬ 
rameters where MI is present, but if and only if 
the MI gain band also contains the zero-frequency 
perturbation as a limiting case (baseband MI). 
We have numerically confirmed that in the baseband-MI 
regime rogue waves can indeed be excited from a noisy in¬ 
put cw background. Otherwise, when there is passband 
MI we only observed the generation of nonlinear wave 
oscillations. Based on the above findings, we are led to 
believe that the conditions for simultaneous rogue wave 
existence and of baseband MI may also be extended to 
other relevant and integrable and non-integrable physi- 
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cal models of great interest for applications, for instance 
consider frequency conversion models [SH [35] where ex¬ 
treme wave events and complex breaking beaviours are 
known to place [36l |37j . 
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